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APPENDIX TO "Target Ligand Generation" 

Naming Inventors: Albert Pierce and Guy Bemis 

1 # ! /usr/bin/python2 . 2 
2 

3 import sys , string , math 

4 from openeye . oechem import * 
5 

6 worstDist = 1.0 #max distance b/w bonds across which molecules will 

7 be crossed 

8 worstAng = 15.0 #max angle b/w bonds across which molecules will be 

9 crossed 

10 worstDist = worstDist**2 
11 

12 def bondmatch (bl, b2) : 

13 #returns "matchsense" = 1 if Moll .bondl . Begin corresponds to 

14 Mol2 . bond2 . Begin 

15- # = 2 if Moll .bondl .Begin corresponds to 

16 Mol2 .bond2 .End 

17 # matchsense = 0 if bonds aren't close (+/-worstDist) and 

18 parallel ( + /-worstAng) 

19 # = 0 if bonds are of different order 

20 global worstDist 
" 21 

22 parallel = 0 

23 antiparallel - 0 

24 if bl.GetOrderO »= b2 .GetOrder ( ) : 

25 return 0 

26 at 11 = bl.GetBgnO 

27 at 12 = bl.GetEndO 

28 at 21 = b2.GetBgn() 

29 at22 = b2.GetEnd() 

30 cartll = moll .GetCoords (atll) 

31 cartl2 = moll . GetCoords (atl2 ) 

32 cart21 = mol2 . GetCoords (at21) 

33 cart22 = mol2 . GetCoords (at22 ) 
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34 

35 distl = r2 (cartll, cart21) 

36 if distl<worstDist : 

37 dist2 = r2 (cartl2, cart22) 

38 if dist2<worstDist : 

39 if deadend (atll) and deadend (at21) : 

40 return 0 

41 if deadend (atl2) and deadend (at22 ) : 

42 return 0 

43 parallel = 1 
44 

45 distl - r2 (cartll, cart22) 

46 if distl<worstDist : 

47 dist2 = r2 (cartl2 , cart21) 

48 if dist2<worstDist : 

49 if deadend(atll) and deadend (at22 ) : 

50 return 0 

51 if deadend (atl2) and deadend (at 21) : 

52 return 0 

53 antiparallel = 1 
54 

55 if (parallel+antiparallel) ==0 : 

56 return 0 
57 

58 ang = calcAngl (cartll , cartl2 , cart21 , cart22 ) 

59 if ang > worstAng: 

60 return 0 

61 if parallel==l: 

62 return 1 

63 else: 

64 return 2 
65 

66 def deadend (atom) : 

67 count = 0 

68 for nbor in atom . Get Atoms () : 

69 count += 1 

70 if count < 2: 

71 return 1 
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72 return 0 
73 

74 def r2(xyzl, xyz2): 

75 dist2 = (xyzl[0] - xyz2[0])**2 

76 dist2 = dist2 + (xyzl[l] - xyz2 [1] ) **2 

77 dist2 = dist2 + (xyzl[2] - xyz2[2])**2 

78 return dist2 
79 

80 def calcAngl (xyzl, xyz2 ,xyz3 ,xyz4) : 

81 vectl = [xyzl [0] -xyz2 [0] , xyzl [1] -xyz2 [1] , xyzl [2] -xyz2 [2] ] 

82 vect2 = [xyz3 [0] -xyz4 [0] , xyz3 [1] -xyz4 [1] , xyz3 [2] -xyz4 [2] ] 

83 vect3 = [vectl [0] -vect2 [0] , vectl [1] -vect2 [1} , vectl [2] -vect2 [2] ] 

84 rl = (vectl [0] **2 + vectl [1] **2 + vectl [2] **2 ) **0 . 5 

85 r2 = (vect2[0]**2 + vect2 [1] **2 + vect2 [2] **2) **0 . 5 

86 r3 = (vect3[0]**2 + vect3 [1] **2 + vect3 [2] **2 ) **0 . 5 

87 angl = math . acos ( (rl**2 + r2**2 - r3**2) / (2*rl*r2) ) 

88 angl = (angl /math. pi) *180 

89 if angl>90: 

90 angl = 180-angl 

91 return angl 
92 

93 def avgCoords (atoml , atom2 ) : 

94 i = 0 

95 cart - [] 

96 coords2 = atom2 . GetParent ( ) .GetCoords (atom2) 

97 for coord in atoml . GetParent ( ) .GetCoords (atoml) : 

98 cart . append ( (coord + coords2 [i] ) /2 ) 

99 i+=l 

100 return cart 
101 

102 def GetNbrs (atom) : 

103 AtmList = [] 

104 for bond in atom .GetBonds ( ) : 

105 AtmList . append (bond. GetNbr (atom) ) 

106 return AtmList 
107 

108 def splitmol (mol , bond, otheratom, otheratom2 ) : 

109 newmoll = OECreateOEMol (mol) #Initiate new molecule 

46 



14435-004001 /VPI 03-12 



APPENDIX 



110 newmol2 = OECreateOEMol (mol) #Initiate new molecule 
111 

112 atoml a newmoll. GetAtom (HasAtomldx (bond. GetBgnldx ()) ) 

113 atom2 = newmoll .GetAtom (HasAtomldx ( bond. Get End I dx ( ) ) ) 

114 atoml .GetParent () .SetCoords (atoml , avgCoords (atoml , otheratom) ) 

115 atoml. SetName ("Al") 
116 

117 atomlist = [atom2 . Getldx ( ) ] 

118 newatoms = atomlist 

119 forbidden = [atoml . Getldx () ] 
120 

v 121 while len (newatoms) : ftrecursively remove neighboring 

122 atoms/bonds from molecule 

123 for atom in newatoms: 

124 newatoms . remove (atom) 

125 Atom = newmoll. Get Atom (HasAtomldx (atom) ) 

126 if Atom is None: continue 
127 

128 AtmList = GetNbrs (Atom) 

129 for nbor in AtmList: 

130 if nbor .Getldx () in forbidden: 

131 Atom2 = newmoll .GetAtom (HasAtomldx (forbidden [0] ) ) 

132 newmoll .Dele teBond (newmoll . GetBond (Atom, Atom2) ) 

133 continue 

134 newatoms . insert ( 0 , nbor . Getldx ( ) ) 
135 

136 for rmbond in Atom. GetBonds ( ) : 

137 newmoll .DeleteBond (rmbond) 
138 

139 newmoll .Dele teAtom (Atom) 
140 

141 mollist = [newmoll] 
142 

143 atoml = newmol2 . GetAtom (HasAtomldx (bond . GetEndldx ()) ) 

144 atom2 = newmol2 . GetAtom (HasAtomldx (bond. GetBgnldx ()) ) 
145 

146 atoml. GetParent () .SetCoords (atoml , avgCoords (atoml , otheratom2 ) ) 

147 atoml. Se tName ("Al") 
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149 atomlist = [atom2 . Getldx ( ) ] 

150 newatoms = atomlist 

151 forbidden = [atoml . Getldx () ] 
152 

153 while len (newatoms ) : #recursively remove neighboring 

154 atoms/bonds from molecule 

155 for atom in newatoms: 

> 156 newatoms . remove (atom) 

157 Atom = newmol2 .GetAtom (HasAtomldx (atom) ) 

158 if Atom is None: continue 
159 

160 AtmList = GetNbrs (Atom) 

161 for nbor in AtmList: 

162 if nbor. Getldx () in forbidden: 

163 Atom2 = newmol2 . GetAtom (HasAtomldx (forbidden [0] ) ) 

164 newmol2 .DeleteBond (newmol2 .GetBond (Atom, Atom2) ) 

165 continue 

166 newatoms . insert ( 0 , nbor . Getldx ( ) ) 
167 

168 for rmbond in At om. Get Bonds ( ) : 

169 newmol2 .DeleteBond (rmbond) 
170 

171 newmol2 . DeleteAtom (Atom) 
172 

173 mollist . append (newmol2) 

174 return mollist 
175 

176 def linkmol (mol , bondorder) : 

177 atomlist = [] 

178 for atom in mol . Get Atoms ( ) : 

179 if atom.GetNameO == "Al": 

180 atomlist . append (atom) 

181 if len (atomlist) == 2: 

182 mol .NewBond (atomlist [0] , atomlist [1] , bondorder) 

183 return mol 
184 

185 if len (sys . argv) ==3 : 
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186 output = sys.argv[l] 

187 if (output != '-osdf') and (output 1= 1 -osmi ' ) : 

188 print "usage : breed. py [-osdf | -osmi] file.sdf" 

189 sys.exit(l) 

190 else: 

191 print"usage: breed. py [-osdf | -osmi] file.sdf" 

192 sys.exit(l) 
193 

194 ifs = oemolistreamO 

195 ofs = oemolostreamO 

196 of s .SetFormat (OEFormat_SDF) 

197 # Read SD file 

198 mollist = [] 

199 smilist = {} 

200 if (ifs.open(sys.argv[2] ) ==1): 

201 for mol in if s .GetOEMols ( ) : 

202 OESuppressHydrogens (mol) 

203 newmol = OECreateOEMol (mol ) 

204 smi = OECreateCanSmiString (mol) 

205 smilist [smi] = 1 

206 mollist . append (newmol) 

207 else: 

208 print "Error accessing SD file" 

209 sys.exit(O) 

210 # Loop over all bond pairs in all molecule 

211 pairs. 

212 i = 0 

213 for moll in mollist: 

214 i-f-1 

215 for j in range (i , len (mollist )) : 

216 mol2 = mollist [j] 

217 for bondl in moll .GetBonds ( ) : 

218 if bondl.IsInRingO : 

219 continue 

220 else: 

221 for bond2 in mol2 .GetBonds ( ) : 

222 if bond2.IsInRing() : 

223 continue 
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224 match_sense = bondmatch (bondl , bond2 ) #is bond a match b/w 

225 2 molecules? 

226 if (match_sense) : 

227 molname = mollist [i- 1] . GetTitle ( ) + "_ M + 

228 mol2.GetTitle() 
229 

230 if match_sense==l : 

231 atoml = bond2 . GetBgn ( ) 

232 atom2 = bond2 . GetEnd ( ) 

233 else: 

234 atoml = bond2 . GetEnd ( ) 

235 atom2 = bond2 . GetBgn ( ) 
236 

237 fragments - [] 

238 fragments. append (splitmol (mollist [i- 

239 1] , bondl, atoml, atom2) ) #split molecule 1 
240 

241 if match_sense==l : 

242 atoml = bondl .GetBgn () 

243 atom2 = bondl . GetEnd ( ) 

244 else: 

245 atoml = bondl . GetEnd ( ) 

246 atom2 = bondl . GetBgn ( ) 
247 

248 fragments . append ( splitmol (mol2 , bond2 , atoml , atom2 ) ) 

249 #split molecule 2 

250 if match_sense == 1: 

251 OEAddMols (fragments [0] [0] , fragments [1] [1] ) #2nd 

252 half mol 2 

253 bondorder = bondl .GetOrder ( ) 

254 newmoll = linkmol (fragments [0] [0] , bondorder ) 

255 OEAddMols (fragments [1] [0] , fragments [0] [1]) #2nd 

256 half mol 1 

257 newmol2 = linkmol (fragments [1] [0] , bondorder) 

258 elif match_sense == 2: 

259 OEAddMols (fragments [0] [0] , fragments [1] [0] ) #lst 

260 half mol 2 

261 bondorder = bond2 . GetOrder ( ) 
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262 newmoll = linkmol (fragments [0] [0] , bondorder) 

263 OEAddMols (fragments [0] [1] , fragments [1] [1] ) #2nd 

264 half mol 2 

265 newmol2 = linkmol ( fragments [0] [1] , bondorder) 
266 

267 smi = OECreateCanSmiString (newmoll) 

268 if not smilist . has_key (smi) : 

269 smilist [smi] = 1 

270 if (output == • -osdf 1 ): 

271 OEWriteMolecule (of s , newmoll) 

272 else: 

273 print smi,molname 
274 

275 smi = OECreateCanSmiString (newmol2 ) 

276 if not smilist .has_key( smi) : 

277 smilist [smi] = 1 

278 if (output == • -osdf'): 

279 OEWriteMolecule (of s , newmol2 ) 

280 else: 

281 print smi,molname 
282 

283 
284 
285 
286 
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